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Abstract: Spacetime in the vicinity of an event horizon can be probed using obser- 

, vations which explore the dynamics of the accretion disc. Many high energy theories 

D , of gravity lead to modifications of the near horizon regime, potentially providing a 

' testing ground for these theories. In this paper, we explore the impact of braneworld 

gravity on this region by formulating a method of deriving the general behaviour 

r\ \ of the as yet unknown braneworld black hole solution. We use simple bounds to 

c^ ' constrain the solution close to the horizon. 
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1. Introduction 

Black holes are a fascinating topic of study, whether it be to explore the geometry of 
strong gravity and its quantum effects, or to understand the astrophysics of massive 
objects in our galaxy. From the point of view of standard physics in four dimensions, 
black holes are well described by the Kerr-Newman family of solutions, solutions 
in standard four-dimensional Einstein-Maxwell gravity. However, the true nature 
of gravity, though well measured from tabletop to solar system scale experiments, 
is surprisingly less well determined at larger or smaller scales. For example, the 
assumptions of dark matter and/or energy are made in order to fit the observed 
Universe with Einstein theory, yet it is quite possible that it is the theory of gravity 
which should be altered to fit these data instead. 

Theoretically, the possibility that gravity might not be fundamentally four- 
dimensional, or indeed Einsteinian, is gathering credence. This is in part an impact 
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of superstring theory, which is consistent in ten dimensions (or M-theory in eleven di- 
mensions) , but also the more phenomenological recent developments of "braneworld" 
scenarios, JI], 0, have had a direct influence on studies of more exotic gravitational 
theories in four-dimensions. Experimentally, the small scale implications of alterna- 
tive theories of gravity have focused on the impact of Kaluza-Klein (KK) modes or 
small black holes in colliders, [^, Q , cosmic ray showers, [Q , or in particle interactions 
in supernovae or nucleosynthesis, [Q. Large scale implications of modified gravity 
have been conducted mostly within the MOND (modified Newtonian gravity) set-up 
[0], although some preliminary studies have investigated braneworld modifications 
to the microwave background |^. Interestingly, while much theoretical work has 
been done on black holes in modified or higher dimensional gravity theories, there 
is very little direct link with experimental data so far. Yet there are observations of 
astrophysical black holes which provide some constraints on strong gravity. Stellar 
remnant black holes, with masses of about 10 times the mass of our Sun, can form as 
a result of very massive star evolution. Most stars form in binary systems, so if the 
companion star is close enough then it can provide a source of material falling onto 
the black hole. This material has the angular momentum of the orbit so it forms a 
ring, but angular momentum transport (via a magnetic dynamo process, Q) spreads 
this ring out into a disc. The magnetic dynamo also dissipates energy, so this means 
that there is emitting material down to the last stable orbit around the black hole. 
For large mass accretion rates this disc emission takes a rather simple form as the 
material is optically thick. The energy is emitted and absorbed many times before 
it escapes the disc, so thermalizes to a blackbody spectrum at a given radius. The 
total disc spectrum is a sum over all radii of these different temperature blackbodies, 
with the maximum temperature (which is of order 1 keV for these stellar mass black 



holes) being set by the emission from close to the inner edge of the disc [^. The 
combination of observed disc luminosity and temperature can constrain the size of 
the inner disc radius, while dynamical studies of the orbit can constrain the black 
hole mass. Together these give an estimate of the innermost stable orbit in terms 
of Schwarzschild radii. Current observations indicate that this is generally of order 
6GM, as expected for a Schwarzschild black hole, though there are a few objects 
with significantly smaller radii, which are interpreted as moderate spin Kerr black 
holes ^. 



Another observational tracer of the gravitational field comes from a corona above 
the disc, which forms a high energy tail to the disc emission. These high energy X- 
rays illuminate the disc and form an iron fluorescence line, produced from material 
with a large rotational velocity in a strong gravitational field. Thus what is seen at 
infinity is affected by both special (doppler shift, length contraction, time dilation) 
and general (gravitational redshift) relativistic effects. Together these transform an 
initially narrow atomic transition into a broad, skewed and reddened line profile. 
Current observations indicate that the observed lines are consistent with material 
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close to the last stable orbit in Schwarzschild or Kerr black holes [12|. These rel- 
ativistic effects also distort the spectrum of the intrinsic (sum of blackbodies) disc 
emission [ pT^ . 

Thus black holes provide a promising test bed for testing our understanding 
of gravity. While stellar, or galactic, black holes do not provide regions of high 
curvature, at least at the classical level, if we accept as fact that black holes radiate, 
then it is likely that close to the event horizon quantum effects become relevant. 
Alternatively, the effect of higher dimensions or additional fields in the gravity sector 
might begin to make their presence felt as the event horizon is approached. In this 
paper we explore the case of braneworld black holes. 

The braneworld paradigm views our universe as a slice of some higher dimen- 
sional spacetime. Unlike the Kaluza-Klein picture of extra dimensions, where we do 
not notice the extra dimensions because they are so small and our physics is 'av- 
eraged' over them, the braneworld picture can have large, even non-compact, extra 
dimensions which are unobservable at low energies since we are confined to the brane. 
Naturally, at not so low energies, the extra dimensions can have experimental con- 
sequences |jl|, ||. Confinement to the brane, while at first sounding counter- intuitive, 
is in fact a common occurrence. The first braneworld scenarios ^^ used topological 



defects to model the braneworld, and zero-modes on the defects to produce confine- 
ment. Of course in string theory, D-branes have 'confined' gauge theories on their 
worldvolumes. 

The braneworld scenario provides a set-up in which we have standard four- 
dimensional physics confined to the brane, but gravity (and possibly a small number 
of other fields) propagating in the bulk. The new phenomenology of braneworld sce- 
narios is then primarily located in the gravitational sector, with a particularly nice 
possible resolution of the hierarchy problem being its primary motivation. Clearly 
however, the astrophysical and cosmological implications of such a scenario have 
a more immediate, and more directly measurable, impact. For these issues, the 
most popular model to explore, and the one which we will be using, has been the 
Randall-Sundrum scenario, 0, which consists of a domain wall universe living in 
five-dimensional anti-de Sitter (adS) spacetime. This model is loosely motivated by 
the Horava-Witten compactification of M-theory jTEI . 

The Randall-Sundrum model has one (or two) domain walls situated as minimal 
submanifolds in adS spacetime. In its canonical form, the metric of the braneworld 
is 

rfs2 = e-2'=N [rft^ _ ^^2] -dz^ (1.1) 

Here, the spacetime is constructed so that there are four- dimensional flat slices 
stacked along the fifth ^-dimension, which have a z-dependent conformal pre-factor 
known as the warp factor. Since this warp factor has a cusp at 2; = 0, this indicates 
the presence of a domain wall - the braneworld which represents an exactly flat 
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Minkowski universe. 

Randall and Sundrum showed in their original papers that although gravity 
was inherently five-dimensional, and the spacetime was strongly warped, as far as 
a four- dimensional braneworld observer was concerned, the Newtonian potential of 
a particle on the brane was indeed the four-dimensional 1/r potential. These early 
results were backed up by more complete analyses confirming that the graviton prop- 
agator did indeed have the correct tensor structure, and that the effect of the KK 
modes was to introduce a 1/r^ correction to the gravitational potential [|T^ . 

Of course, of real interest in astrophysics and cosmology is not the perturbative 
study of the graviton propagator, but a concrete understanding of the true non- 
perturbative nature of gravity on the brane. This is what is relevant for a black 
hole. To some extent, this was provided by the brane-based approach of Shiromizu 
et. al. (SMS) |]T^, who used a Gauss-Codazzi approach to obtain the braneworld 



gravitational equations: 

G^u = ^4.9 fiu + SttGatT^,^ + K S^y + £^j,y (1.2) 

Here, A4 is a residual cosmological constant on the brane. It represents the mismatch 
between the brane tension and the bulk negative cosmological constant. (In the RS 
scenario, the brane tension was precisely tuned to balance out the negative bulk 
cosmological constant.) The final two terms are the brane gravity effects. The first, 
S^v, consists of squares of the energy- momentum tensor, which is only relevant at 
high energies. The final term, £^^, is a so-called Weyl term, and consists of the 
projection of the bulk Weyl tensor onto the brane. Although the SMS equations 
do provide an apparent nonperturbative description of gravity on the brane, it is 
important to emphasize that the Weyl term is not given in terms of data on the 
brane, and is a complete unknown from the brane point of view. Therefore, in order 
to take a pure brane-based approach to gravity, some assumptions must be made 
about £f^u- 

Two main nonperturbative gravitational problems are of clear interest: Cosmol- 
ogy, and Black Holes. The first problem, that of finding the braneworld generaliza- 
tion of the FRW universe, has been well explored and understood. Brane-based work 
concentrates mostly on an explanation of the "non-conventional" energy momentum 
squared terms, as well as the unknown dark radiation effect from the Weyl term [|1^]. 



For cosmology however, the high degree of symmetry present renders the full five- 
dimensional problem fully integrable |l^, and the general cosmological braneworld 



is fully understood in terms of a slice of a five- dimensional adS black hole ||2^. The 
mass of this bulk black hole then generates the radiation-style Weyl term. It is there- 
fore fair to say that while all the implications may not have been calculated, brane 
cosmology for the pure RS scenario is pretty well understood. 

The situation for the black hole on the brane is somewhat different however. 
Though it would appear that the two cases are similar in that cosmological branes 
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have a dependence on time as well as the bulk ^-coordinate, and black holes on radius 
and bulk 2;-coordinate, in fact, the symmetry groups of the spacetimes are crucially 
different. For cosmology, the metric splits into two parts - the two dimensions on 
which it depends, and the spatial part of the universe, which has constant curvature. 
The mathematics of the cosmological braneworld is therefore a two dimensional field 
theory which turns out to be totally integrable. For the black hole however, the 
metric splits into three parts - the two dimensions on which it depends, the time 
coordinate and the remaining spatial part in which the horizon resides. Thus there 
are two fields in our two dimensional theory, one of which acquires a potential, and 
there is no longer a simple solution |^ . 



The case of the black hole in this braneworld picture now becomes of great 
interest and importance. Of interest as a problem in higher dimensional gravity. 



and of importance because the adS/CFT correspondence, [Q, relates the classical 
five-dimensional braneworld black hole solution to the four-dimensional quantum 
radiating black hole|^, ^. The first attempt, [^, to find a black hole solution 



replaced the Minkowski metric in ( |1 . 1| ) by the Schwarzschild metric, thus creating a 
black string sticking out of the brane. Unfortunately, as suspected by the authors. 



this string is unstable to classical linear perturbations |26|. Chamblin et. al. realised 



that the true localised black hole would be a slice of a five- dimensional accelerating 



black hole metric (known as the C- metric, p^, in four dimensions), however no such 



metric has as yet been found. A lower dimensional version of a black hole living on 
a 2 + 1-dimensional braneworld was however presented by Emparan, Horowitz and 
Myers |^, using this four- dimensional C-metric. Since then, several authors have 



attempted to find the full metric - notably numerical work by Wiseman |^U[, and 
others [^ . 



In this paper, we are interested in near horizon modifications to General Relativ- 
ity, which in the absence of a full five-dimensional solution might seem problematic. 
However, our aim here is not to attempt to answer the full five-dimensional problem, 
but to revisit the four- dimensional braneworld and to take a practical approach to 
finding the braneworld metric. The motivation will be to see if we can categorize 
(preferably analytically) various classes of near horizon behaviour, and to see if there 
is any universality to the four- dimensional braneworld solutions. 

The literature has several special solutions, notably the tidal Reissner- Nordstrom 
solution of Dadhich et. al. |^, but also the PPN parametrized solutions of Casadio 



et. al. [^. Visser and Wiltshire [^] presented a more general method which gener- 
ated an exact solution for a given radial metric form. All of these approaches however 
have in common the assumption of an "area gauge" for the radial coordinate, r, in 
the metric, i.e., that the area of the 2-spheres surrounding the black hole behaves as 
A{r) = Anr"^. Thus the area of spheres surrounding the black hole increases mono- 
tonically between the horizon and spacelike infinity. However, the monotonicity of 
A{r) is only guaranteed if the dominant energy condition holds, and there is no rea- 
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son to suppose that this will be the case for the Weyl term, indeed, Dadhich et. al. 
have Soo < - a violation of the weak energy condition! In fact, energy conditions 
are persistently violated in dimensionally reduced theories of gravity. 

Our reason for suspecting that A{r) is not monotonic lies in the putative higher 
dimensional C-metric, which would consist of an accelerating black hole being 'pulled' 
by a string. The appropriate higher dimensional metric for a Poincare invariant string 
has a turning point in the area function, and the 'horizon' is in fact singular [^. It 
is therefore possible that this renders the black hole horizon also singular. Moreover, 
using the dual description of a CFT with cutoff living on the brane [Q, a static 
quantum black hole must have a singular horizon [^. Since the singularity of the 



string has a diverging area function, it seems likely that the black hole itself might. 
Such a turning point in the area function can be thought of as a wormhole in the 
geometry, so called because the spatial part of the Schwarzschild metric: 



y 'T' I iy^ ry ^ 



(where r = (f + G'M/2)^/f) has this property of the area of spheres decreasing as we 
approach r = 2GM, then stationary at r = 2GM, then increasing again as we move 
onto the other asymptotically flat regime of the maximally extended Schwarzschild 
solution. As we will see, this is a very good analogy, as there is an exact analytic 
braneworld metric which has this spatial form, and simply moves the event horizon 
relative to Schwarzschild either outside r = 2GM or through the wormhole neck onto 
the other Kruskal branch. 

In identifying general features of braneworld solutions, the questions we explore 
in this paper are: 

• When is the horizon singular? 

• When does the area function have a turning point? 

• When are black hole solutions asymptotically flat? 

Within the context of an equation of state of the Weyl term we are able to give 
definitive answers to all of these questions. 

The layout of the paper is as follows. In the next section we derive the braneworld 
equations in a general gauge and show how to reduce these to a two-dimensional 
dynamical system for a given equation of state. At this point we give the analytic 
solution corresponding to the Schwarzschild wormhole. In section ^ we analyze this 
system in general, showing how to answer the questions above. In section ^ we make 
contact with the asymptotic linearized propagator, and work on small black holes, 
and present arguments leading to an analytic near horizon metric. Finally, in section 
^, we use this metric to explore the phenomenology of astrophysical black holes. 
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2. Spherically symmetric braneworld metrics 

We start by looking at the general static, spherically symmetric metric on the 
braneworld. By taking the metric to be static, we are taking the point of view that 
there exists a five- dimensional solution analogous to the C-metric in four- dimensions 
which has a timelike Killing vector, and can therefore be 'sliced' by the braneworld in 
such a way as to create a static four- dimensional black hole on the brane. In this we 
assume that the three-dimensional braneworld black hole constructed from the four- 
dimensional C-metric has a direct analog in one dimension higher. We note that this 
produces a static braneworld black hole, and precludes the possibility of a dynamical 
back-reaction to Hawking radiation, instead corresponding to a Boulware choice of 
boundary conditions, [E3], which give a singular 'horizon'. It should be noted that 



there is not a consensus as to whether the braneworld black hole should be static 
(and therefore singular, by the reasoning of |^), or time dependent, corresponding 
to an evaporating black hole, as explored by Tanaka |]35| . 

2.1 The metric 

The general static spherically symmetric metric on the brane can be written as: 

ds^ = A^{r)df - B\r)dr^ - C\r)dnl (2.1) 

Clearly this is not in the simplest gauge, as we can still choose our radial function, 
r, quite arbitrarily, however, it proves to be convenient to use this over-general form 
in order to choose the best gauge for problem solving, and to compare this to more 
familiar gauges more readily. The main reason for using this form of the metric, with 
an arbitrary function for the area of the 2-spheres rather than r^, is that there is good 
reason to believe that the area function might not be monotonic. With this form 
of the metric, the second derivative of the area radius, C, {i.e., the radial function 
defined by y^A/^vr) is given by the following combination of the Einstein tensor: 

therefore, for the area function to be guaranteed to be monotonic, we must have 
Gl — C!^ > 0, which is equivalent to the dominant energy condition. In usual Einstein 
gravity this is generally the case, but when there are extra dimensions, or extra fields. 



this is no longer the case. For example, in the cosmic p-brane, |3J] the area function 
actually blows up on the 'horizon'. Since the cosmic p-brane might be expected to 
have some connection to the higher dimensional C-metric, possibly rendering the 
'horizon' singular, it is vital that in any exploration of braneworld black holes we do 
not make the restrictive ansatz of C = r. 

The 'vacuum' brane equations from (|1.2| ) are 
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We follow Maartens [3^ in using the symmetry of the physical set-up to put the 
Weyl energy into the form: 



£^^ = U {u^Uu - Ih^^) + n (r^r^ + ^h^^) 



[2A) 



where u^ is a unit time vector, and r'^ a unit radial vector. Note that the Weyl 
energy is in 'Planck' units, i.e., there is no preceding SirG since S is derived from 
gravitation in the bulk. If we want to compare with 4D matter, we should rescale 
by I/SttG. 

We now have the equations of motion: 

{u + 2n) 



Gl 



1 



1 
52 



C" C'B' 

2 2 

C CB 



U 



r^ = J- - ^ 



AC C"2 



(-it) _ 



A" C" 



A'C A'B' B'C 



AC 



AB 



BC 



(w-n) 



(2.5) 
(2.6) 

(2.7) 



An alternate and useful equation is the Bianchi identity: 

{u + 2n)' + 2^ (n + 2U) + en^ 







(2. 



(or conservation of "energy-momentum"). 

This system of equations has been solved in many special cases, and more general 
techniques have also been presented. Briefly, the special cases are the tidal Reissner- 
Nordstrom solution of Dadhich et. al. |HT|, or the solutions which assume a given form 
for the time or radial part of the metric [^. Visser and Wiltshire [Q presented 
a more general method which generated an exact solution for a given radial metric 
form. In all of these cases however, the radial gauge C = r was chosen (although 
[0 did comment on how to use their method when C was not monotonic). 

A reasonable alternative to making guesses for various of the metric functions is 
instead to follow a pragmatic approach as in cosmology. When solving for an FRW 
universe, the precise details of the composition of the universe are approximated 
by an isotropic perfect fluid energy-momentum tensor, and, most pertinently, an 
equation of state is assumed for this source. Dust {p = 0) for the later universe, 
radiation {p = p/3) for somewhat earlier times, and cosmological constant {p = —p) 
for inflation. Clearly the actual evolution and matter content of the universe is 
more detailed and complicated, but this method is useful, accurate, and universally 
accepted. Here we propose the analogous approach: an equation of state for U and 
n, n = ^^U. Of course, a priori there is no reason to suppose that the Weyl 
term should obey an equation of state (and indeed we will argue later that it will 
not) however, it is quite possible that just as the universe is approximately matter 
dominated or radiation dominated in certain eras, the Weyl term might have certain 



asymptotic equations of state which may be useful as near-horizon or long range 
approximations. 

2.2 The dynamical system 

In the next section we provide a detailed analysis of all equations of state, using the 
common method of dynamical systems analysis of the equations of motion. In our 
case, it is clear that the most convenient gauge for this type of analysis is B = C. 
In this case, writing 

n' A' r" 



we have: 



2X' + X^ - 1 = 3UC^ = -{XY - 1) (2.10) 

7 

1-2Y' -Y^ =UC^ = -{XY -1) (2.11) 

7 

together with the constraint 

l-XY = -^{U + 2n) = -7^ (2.12) 

The plot of the {X, Y} phase plane gives us curves which are solutions, X{p), Y{p), 
to this dynamical system. Whether or not these trajectories correspond to an actual 
black hole depends on whether the integrated solutions A{p), C{p) have the behaviour 
we expect for a black hole, such as an horizon, asymptotic flatness and so on. 

It is therefore useful before proceeding with the analysis to extract some asymp- 
totic information from the Einstein equations in order to identify general features of 
the solution on the phase plane. 

2.3 Asymptotic analysis 

There are two clear asymptotic regions in which we would like to have some infor- 
mation about the black hole solution - the horizon and infinity. Near the horizon 
A ^ 0, and we therefore expect that X will become infinite. Near infinity however 
we would like spacetime to be asymptotically flat, which means (in the area gauge) 
A ~ 1 - 0{r-^), 5 ~ 1 - 0{r-^), C = r. Or A ~ 1 - 0(6"^), 5 = C ~ e'' + 0(1) in 
dynamical system (d.s.) gauge. 

For the purpose of clarity, we will look at the far field region in the area gauge, 
in which the leading behaviour of the Weyl energy and radial metric component is 
easily read off as: 

U = ^ (2.13) 

5- = 1 - ^ - ^ (2.14) 

r 3r 
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What this clearly shows is that for < 7 < 3 the braneworld metric does not 
tend to Minkowski spacetime for large r. Indeed, even for 7 > 3 the metric is 
not asymptotically flat, since asymptotic flatness requires corrections to Minkowski 
spacetime to be 0(l/r). Therefore, an asymptotically flat braneworld black hole 
requires an asymptotic equation of state 7 < 0. 

Near the horizon on the other hand, we cannot assume the area gauge, and 
integrating the energy conservation equation gives 



3 



U oc A-^r^C"^ (2.15) 

Whether or not the horizon corresponds to a singularity in the Weyl energy depends 
on the magnitude and sign of 7, as well as on the behaviour of C. 

2.4 A simple analytic solution 

We conclude this section by demonstrating the use of both the general gauge C{r), 
as well as the dynamical system by deriving an alternate form of a known analytic 
solution with U = 0. A quick glance at the equation of state shows that this is 
formally the limit 7 — > ±00. Since letting 7 become infinite decouples the equation 
of state from the dynamical system, the phase plane is in fact the same whether 7 is 
+00 or —00. Our main reason for deriving this solution afresh is that we wish to use 
it in section ^ for a near horizon limit, but clearly, any exact solution is helpful, and 
our derivation makes the overall structure of the spacetime somewhat clearer than 
using an area gauge. 

For U = 0, (|2lO| , ^ll| ) have the solutions: 



X,Y = tanh ^-P^ , coth ^^ (2.16) 

Clearly the former solution is appropriate for Y, as this corresponds to 

C = Cocosh'^-^—^ (2.17) 

which, since we expect an asymptotically flat spacetime which has r = e^, gives 
Co = 4e^o . We then have for A: 

21nA + ln(g = 21ncosh ^^~^^^ (2.18) 



or 



I.e. 



21nA + ln(:7 = 21nsinh^^^ — ^ (2.19) 



A-A,-C, cosh^ ^^^^ = cosh^ ^^^^ or sinh^ ^^^^ (2.20) 
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For an horizon, we expect that A^ will have a zero, hence we choose the second 
solution, giving 



^2 ^ ^ /cosh(p - pi) - 1 



Co \cosh(p-po) + 1 



(2.21) 



For asymptotic flatness Aq = 4e^^ 
Overall the metric is: 



Vcosh(p-po) + iy 2 L^ II} K ) 



or, writing r = e'': 



ds' 



(r + ro)^ 
this metric has appeared in the area gauge as 



Z2 ""-fJ 



:de- 



-dr'- 



-dVti 



(2.23) 



ds' 



, ,^ 2GM 



df 



dH 



1- 



2GM ^ ^^11 
R 



(2.24) 



where i? = (r + ro)^/r, GM = 2ro, and GMe = ri — tq. The anisotropic stress for 
this solution is 11 = 3GMe/AC^. 





Figure 1: Pictorial representation of a constant time slice of the metric ( |2.23| ). 

While the area gauge gives a familiar spatial part of the metric, note that at 
the old Schwarzschild 'horizon', R = 2GM, the t — t component of the metric does 
not vanish. For e > 0, gu will be zero before r = 2GM, and the area gauge holds 
outside the black hole. (Strictly of course, this is not a black hole, as the 'horizon' 
is singular.) If e < then gu = at r < 2GM, and r = 2GM becomes a coordinate 
singularity accessible by timelike observers, the result of choosing an inappropriate 
gauge. 
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3. The general dynamical system 

In order to find a general solution for the braneworld black hole, subject to the 
proviso that the Weyl energy obeys an equation of state, we will now analyze the 
general dynamical system: 

1 -X^ 9 

x' = ^— + ^(xr-i) (3.1) 

1 -y2 Q 
y' = ^— + -(i-xr) (3.2) 

Once again we emphasize that we do not believe that the equation of state will 
necessarily hold over the whole horizon exterior, simply that it will potentially give 
various asymptotic behaviours near infinity or near the horizon. 

The plot of the {X, Y} phase plane gives us curves which are solutions, X{p), Y{p), 
to the dynamical system (|3.1| , p.2|) . These in turn can be integrated to find A{p), C{p). 
A selection of representative phase plane plots are given in figure ||. What these plots 
show is that there are attractors in the phase plane, and that a general trajectory 
flows in from inflnity to one of these attractors. Whether or not these trajecto- 
ries correspond to an actual black hole depends on whether the integrated solutions 
A[p), C{p) have the behaviour we expect for a black hole, such as an horizon, asymp- 
totic flatness and so on. 

Briefly, a horizon corresponds to A —>■ 0, therefore we expect X becomes inflnite 
for the black hole horizon. Flat space on the other hand corresponds to A' = 0, 
C oc e^, hence X = Y = 1. A black hole solution must therefore be a trajectory 
which comes in from large X (possibly large Y as well) and terminates on (1, 1). To 
flnd out which equations of state allow this, and whether the area function has any 
turning points, which correspond to zeros of Y, requires a detailed analysis of the 
phase plane. 

3.1 Features of the phase plane 

Part of any dynamical systems analysis is an identiflcation and classification of critical 
points, invariant submanifolds, and other distinguishing features of the phase plane. 
In this case, the invariant hyperboloid is easily identified from the constraint ( p.l2| ) 
as XY = 1. Along this hyperboloid {U + 211) = 0, therefore, apart from the special 
case 7 = 0, this corresponds to a vanishing Weyl term, and hence the pure Einstein 
case - i.e., the Schwarzschild solution. This is given by ( p.l7] , p.2lD with po = Pi and 
v4q = Cq. If 7 = 0, in addition to the Schwarzschild solution we have X = ±v3 = 
1/y, which is covered by the analysis of the critical points below. 

Critical points 

The system has 4 critical points {X' = Y' = 0): 

P± = ±(l,l) (3.3) 



-12- 




y = 7 






Figure 2: Phase plane for a range of equations of state. 



Q^ = ±{2i+^^r^'\9-^,^+^) 



(3.4) 



The Q± critical points move as 7 runs from +00 to —00, from Q± = (=1=1, ±1) to 
Q± = (±1, =f1)- For 7 = 3, P and Q are coincident. The nature of the critical points 
is as follows: 

• 3 < 7 P+ is an attractor and P_ a repellor, the Q's are saddle points. 

• < 7 < 3 (5+ is an attractor and Q- a repellor, the P's are saddle points. 

• 7 < P+ is an attractor and P_ a repellor, the Q's are saddle points. 
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The critical point (1,1) corresponds to flat space: 

X = Y = 1^A' = 0, C = C 

^C = eP (3.5) 

Therefore any asymptotically flat solution must terminate on this critical point. 
While P+ is an attractor, this is not really a problem, but for the range of 7 where it 
is a saddle point, only the invariant hyperboloid can satisfy this, and by deflnition, 
this is where we have the exact Schwarzschild solution. 

The critical point (5+ on the other hand, corresponds to a non-asymptotically 
flat spacetime, which, for 7 > —3 in area gauge is: 

(3 + 7) 

which we do not expect to be appropriate to the metric for an isolated source. This 
solution can also be used for 7 < —3, provided one remembers that increasing r 
actually corresponds to moving towards the black hole {gu — ^^ as r -^ 00). This is a 
genuine wormhole solution, in that the area A{r) increases unboundedly towards the 
event horizon, which is located at inflnite proper distance. Unlike the Schwarzschild 
wormhole however, this inner asymptotic region is not flat, but, as already mentioned, 
leads in to a null asymptopia. It is perhaps worth noting that the critical value 
7 = —3 corresponds to the marginal case of no wormhole, but an inflnite throat: 

ds^ = r^dt'-^-rldnl (3.7) 

exactly analogous to the extreme Reissner-Nordstrom black hole throat. 

The general solution of the braneworld black hole therefore requires 7 > 3 or 
7 < in order to terminate on the critical point P+. Luckily, this range of 7 is 
precisely that for which spacetime asymptotes flat space, although only for 7 < is 
it actually asymptotically flat. 

Asymptotes 

Finally, the other region of interest, the black hole horizon, corresponds in general 
to large values of X, for which we can identify the characteristic behaviour. 

The line F = (7 — 3)X/(7 + 9) is an asymptote for 7 G [0, 3], but a separatrix 
for negative 7, and 7 > 3. It corresponds to the solution 

247 47(7-3) 

ds^ = {p- po)^^^^dt^ - (p - po)^^+^ [dp^ + dnl] (3.8) 

which has a singular horizon. 

For negative 7, and 7 > 3 the asymptotic solution is X = 2/(p — po), and 
Y oc {p — Pq)~^''. This gives the metric 

ds' = (p - p,fde - Cl (1 +c,{p- p,f-'/^) [dp' + dn]j\ (3.9) 

The horizon is singular in this case for I7I > 3. 
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3.2 Special solutions 

Because of the number and nature of the critical points, there can be special solu- 
tions which start on one critical point and terminate on another. Depending on the 
value of 7, the attractor P+ can have up to three special solutions corresponding to 
trajectories from each of the other critical points. It is easiest to demonstrate these 
solutions for the extreme case 7 = cxo, where we have an analytic solution of the 
phase plane. 

• P_ solution: 

The solution from P_ to P+ is easy to define in the analytic case: it is X = 
tanh(p — pi)/2, Y = tanh(p — po)/2, which corresponds to a two parameter family 
solution of metrics: 

ds^ = j^^c^t^ - ^^^dr^ - ^^^dnl (3.10) 

(r + To) r* r^ 

For ro = Ti, this is the Schwarzschild wormhole, and corresponds to the straight line 
path between P_ and P+. For ri ^ Tq, the solution corresponds to the other paths 
shown in figure ^ between the critical points. The wormhole is supported here by the 
negative anisotropic stress 11 = —6/C^ = — 6r^/(r + tq)^. See |Q for more detailed 
discussions of wormhole solutions on braneworlds. Although for finite 7 the form 
of this metric will change, the general nature - that of a wormhole connecting two 
asymptotically fiat regions - will not. This solution exists for 7 > 9 and 7 < —3. 

• Q- solution: 

This solution only exists in the 7 = cx) limit. It has X = 1 and Y = tanhp/2 
which integrates to the metric 

ds^ = , '' df - ^'^'/ dr^ - ^'^'/ dnl (3.11) 



(r + tqY ^ 



which can be seen to be a limiting form of the metric ( p.23| ) i.e., a wormhole with 



an inner asymptotic null asymptopia. For 7 < —3, the Q+ critical point has F < 
and hence the (5+ -^ P+ solution actually takes this qualitative form. 

• (5+ solution: 

This solution has F = 1 and X = tanh p/2 for the analytic form. This integrates 
to the metric 

ds^ = (l + -\ dt" - dr^ - r^dQJj (3.12) 

which is not a wormhole, but a flat 3-space with a distorted Newtonian potential. 
The solution will have this form for 7 > 3, but for 7 < —3, the Q^ critical point 
moves below the X-axis, and the solution joining (5+ to P^ takes the above form 
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( |3.11 ). For 7 = —3 this trajectory is identical to the extremal Reissner- Nordstrom 



solution: 



^0\ J-^2 /i ^0\ ,2 „2jn2 



ds' = (1 - -j dt' - (1 - -j dr' - r'dnjj (3.13) 

In addition, there are special values of 7 for which there are Q- —>■ Q+ trajecto- 
ries. Specifically, 7 = 9 has a wormhole: 

a 3/r + roV,2 (^ + ^o)'^2 



ds' = ^^^^.dt^ - i ^^J dr^ - '-^^d^n (3.14) 

and 7 = —3 has a 'throat' solution with bouncing Newtonian potential: 

ds^ = Al cosh^ r df - dr^ - Rldnjj (3.15) 

or adS2 x S"^. 

3.3 Summary 

To summarize: the phase plane shows clearly that flat spacetime is a critical point, 
which is an attractor for equations of state with 7 > 3 or 7 < 0. Only for this 
latter range of 7 is the spacetime asymptotically flat. For these ranges of 7 the 
near horizon behaviour is also given by a defined class of metrics ( |3.9| ), which are 
singular for I7I > 3. Solutions are allowed both with and without turning points in 
the area function, though only for 7 < —3 do these correspond to solutions which 
are asymptotically fiat. 

In addition, there are special solutions which correspond to trajectories between 
critical points. The trajectory from P_ to P+ represents a solution with two asymp- 
totic fiat regions, and is a genuine braneworld wormhole. The solution from one Q 
critical point also represents a wormhole spacetime, which in the case of 7 = —3, has 
a throat. 

4. Physical solutions 

In order to decide what a reasonable black hole metric might look like, it is useful 
to compare to known results, that of the linearized metric, and also the small black 
hole approximation. 

4.1 Comparison with linearized propagator 

One region in which we do know the metric of the black hole is at large r. Randall 
and Sundrum, [^ , computed the corrections to the Newtonian potential, and Garriga 
and Tanaka, [0, the correct tensor propagator: 

2GM ( 2/2\ 1GM ( 1/2 



^« = l-^(l+3-2) ' ^^. = -1-^(1 + 3-2) (4-1) 



-16- 



where we write r to distinguish from the r-coordinate of the area gauge. Transforming 
to area gauge gives: 

2GM AGMP -i_r-2_ 2GM 2GMl'^ 



gu = A^ = l -^— ; -g;; = B-' = l — (4.2) 



Substituting into the equations of motion gives an equation of state 11 = 5W/4, or 
7 = -3/2. 

Interestingly, this asymptotically fiat solution holds down to the event horizon, 
giving a fully nonsingular solution. The Weyl energy: 

-AGMP ( 2GM 4GMl^\ , , 

has a maximum at r ~ 2.4G'M and returns to zero at the event horizon. This is not 
the behaviour we would expect from the Weyl energy, hence 7 = —3/2 should be 
taken only as an asymptotic equation of state. 

4.2 Comparison with small black holes 

Another common limit used in braneworld black holes is the small black hole limit. 
Here, if the black hole has mass much less than the adS scale, it is expected that the 
adS curvature has very little effect on the event horizon, and so the five-dimensional 
Schwarzschild solution in a good approximation, i.e. B~^ = A^ = 1 — -^ in area gauge. 
A quick check of the asymptotic solution (|2.14| ) shows that 7 = —3, or 11 = — 2W. 
This has U = —fi/r^. 

A plot of the phase plane for 7 = — 3 shows that these trajectories form a 
stable family of solutions terminating on P+. There are no trajectories crossing the 
X-axis as y = is a solution to the equations of motion. In other words, there are 
solutions to the black hole system which represent truly infinite throats with varying 
Newtonian potentials. For 7 < —3 once again there are trajectories which cross the 
F-axis, and hence have a bounce in the area function, or a wormhole. 

4.3 A possible black hole solution 

Having compared the equation of state method with known asymptotic solutions, it 
is clear that if an equation of state applies, it is generally a negative one. Moreover, 
from the intuition gained by looking at very small black holes, 7 curiously becomes 
more negative closer to the horizon. What sort of equation of state might hold near 
the event horizon? To explore this, we return to the holographic correspondence. 

In string theory, it has been realized for some time that there is a correspon- 
dence between string theory on adS space, and a CFT on the boundary of that adS 



space |22[. In other words, all of the information contained in the five-dimensional 
gravitational spacetime is encoded in a pure quantum field theory (no gravity) living 
on a four- dimensional spacetime. In the braneworld picture, the brane is not at the 
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Figure 3: Phase plane for 7 



adS boundary, but at a finite distance, and the theory on this brane now contains 
gravity, as well as a conformal energy- momentum tensor - the Weyl term £^y. The 
effect of the brane on the adS/CFT correspondence therefore is that the bulk theory 
of gravity in five dimensions corresponds to the four dimensional brane theory of 
a CFT with a UV cutoff interacting with gravity pSl. Since the brane theory is 



a quantum theory, the holographic correspondence indicates that the classical bulk 
solution projects to a quantum corrected solution on the brane [p3| . 

For cosmological solutions, there is a nice holographic interpretation, where the 
braneworld cosmology can have a radiation source which is the result of the projection 
of the Weyl curvature of the black hole on the brane. This source can be interpreted 
as a CFT in a thermal state corresponding to the Hawking temperature of the bulk 
black hole. The brane cosmological metric has a constant curvature spatial part, and 
its symmetries demand that only a radiation energy Weyl term is allowed. From the 
bulk perspective, this means that every point on the brane is at the same distance 
from the bulk black hole. Thus a flat universe corresponds to a 'fiat' bulk black hole, 
a closed universe to a conventional spherical bulk black hole. One way of imagining a 
braneworld black hole forming is to transport this bulk black hole in towards and onto 
the brane. In doing this, one breaks this symmetry of equidistance from the black 
hole. In other words, from the holographically dual brane point of view, we introduce 
an anisotropic stress by shifting the black hole from its equidistant position. The 
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closer we bring the black hole towards the brane, the more anisotropic the set-up. 
Therefore, we might expect that 11 becomes more and more important both as we 
transport the black hole towards the brane, or as we move closer to the event horizon. 
Therefore a physically reasonable expectation might be that equations of state with 
large 7 are relevant near the horizon. 

Taking this reasoning to its logical extreme, we therefore propose as a "working 
metric" for the near horizon solution the analytic form of section |2.4]: 



d.= ^|^*=-i^*.^-i^.«?, (4.4) 

Clearly, although there is a degree of arbitrariness in this choice, we believe 
that the horizon is likely to be singular, and also that a turning point in the area 
function is also likely. This metric has both these features and the added advantage 
of analyticity, which means that many properties can be calculated explicitly. 

5. Testing the solution and future horizons 



We now wish to explore some simple tests of the near horizon form (|2.23|/^I^) or 



ds' 






2 



^2 dR D2jn2 



dt--. 2GM-Rdnjj (5.1) 



R 

Although this latter form of the metric is not valid throughout the exterior of 
the horizon if e < 0, it is an easier gauge to see the contrasts with the Schwarzschild 
metric, and since it turns out that there can be no stable orbit within the wormhole 
region, it is quite satisfactory for working with accretion discs. 

The first point to note about (|5.1|) is that the ADM and gravitational mass 
(defined by Qu) are no longer the same. Once again, this is quite a common occurance 
in string gravity, and is a result of the extra degrees of freedom of the gravitational 
field. As a result, the weak field tests of light bending and perihelion precession will 
be modified at the 0(e) level, in a way consistent with the PPN formalism. 

The near horizon metric (2.24 in area gauge or 4.4 in isotropic coordinates) is very 
different to the standard Schwarzschild solution. The horizon is always singular even 
for apparently negligibly small e (except for e = which corresponds to S^u = /, so is 
identically Schwarzschild). This is a true singularity, as the energy density from S/^iu 
becomes infinite at the horizon. An intrepid observer plunging into a supermassive 
black hole, expecting to sail seamlessly through the horizon and explore spacetime 
close to the singularity before finally succumbing to tidal forces, is instead crushed 
out of existence at the horizon. Indeed, they are slammed into the horizon at infinite 
proper speed! The are no timelike or null geodesies which connect anywhere outside 
the horizon with the standard Schwarzschild singularity at r = 0. Infalling matter 
or light simply cannot reach this point. 
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e -0.5 -0.1 0.1 0.5 1 2 

~^h oo* 2.02* 2 2m 2^25 2^67 3.6 

Rms 4.26 5.63 6 6.37 7.88 9.82 13.75 

h^ 4.15 10.11 12 14.03 23.78 39.54 83.06 

Rph 2.25 2.83 3 3.17 3.91 4.86 6.82 

Table 1: Radii for the horizon, Rh, minimum stable orbit for particles, Rms and the 
unstable photon orbit Rph in units of GM. We also give the associated angular momentum 
of the minimum stable orbit, /i^. Note that for negative e, the star indicates that the 
horizon is on the lower Kruskal branch, and therefore inside R=2GM. 



Such a drastic change to the spacetime surely has observable effects above the 
horizon. The metric (4.4) can easily be compared with the PPN formalism, and 
observational limits on these parameters from solar system tests set limits on |e| < 
10~^. However, these derive from solar system tests in the weak gravitational field 
regime. This would only be applicable if the metric (4.4) covered the entire spacetime, 
but as discussed above we envisage this only as the near horizon asympote of a more 
general metric. The only mathematical limit on the near horizon metric is e > — | 
from the requirement that ri is positive definite. 

Here instead we constrain e by connection to the observations of accreting black 
holes in our galaxy. The luminosity and temperature of the blackbody emission from 
the accretion disc can be measured from X-ray spectra, giving a direct estimate of 
the emitting area. There are uncertainties in this approach, but some confidence 
can be derived from the fact that changes in the luminosity give rise to the L cc T^ 
behaviour expected from a constant emitting area. Assuming that the uncertainties 
are less than a factor of 2 then the data strongly constrain the minimum stable orbit 
to be < 12G'M. 

We assume that the 'near horizon' metric applies on these size scales and solve the 
Euler-Lagrange equations from (2.24) in the equatorial plane to find the innermost 
stable particle orbit, and its associated angular momentum h = \g,j)^\(f). Table 1 
shows these as a function of e, together with the horizon position {gu = 0) and the 
radius at which light can (unstably) orbit around the black hole. We can use the 
area gauge even for cases where e < (which have a wormhole) as the coordinate 
singularity aX R = 2GM is below the radii of interest for light and particle orbits. 
The observations of accretion disc size in X-ray binary systems give an upper limit 
to e of about 2, and are easily consistent with our expectation of e < 0. (In the 
absence of a complete metric from larger to smaller r, we have taken M to be the 
ADM mass.) 

If this asymptotic metric holds to within a few tens of Schwarzschild radii of the 
horizon then there are potentially observable effects on the spacetime around stellar 
mass black holes. We will explore these further in a subsequent paper. 
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Although we have focused on braneworlds, the techniques will clearly be applica- 
ble to stringy black holes, or indeed any theory that modifies gravity. In the absence 
of a concrete solution to test, our approach of trying to find classes of behaviour 
and some sort of universality near the event horizon seems the best way to proceed. 
Typically, any covariant theory which modifies the near horizon spacetime structure 
of the black hole will have these qualitative features, hence similar effects that we 
have been exploring. All in all, this seems a fruitful alternative to explore in terms 
of testing ideas in high energy gravity. 
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